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ABSTRACT 

We investigate unitary one-matrix models coupled to bosonic quarks. We de- 
rive a flow equation for the square-root of the specific heat as a function of the 
renormalized quark mass. We show numerically that the flows have a finite number 
of solitary waves, and we postulate that their number equals the number of quark 
flavors. We also study the nonperturbative behavior of this theory and show that 
as the number of flavors diverges, the flow does not reach two-dimensional gravity. 
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1. Introduction 



The study of matrix models has led to the uncovering of many rich structures. 
Not only do they describe two dimensional gravity coupled to conformal matter [1 31 , 
at least perturbatively, but they have also been shown to satisfy a KdV hierarchy 
for the Hermitian case' 4 ' 51 , and a mKdV hierarchy for the unitary case' 6 ' 71 . 

The generic one-matrix model satisfies the differential equation' ' 81 



x = J2(k + l/2)T k R k [v], (1.1) 

where the Rk[u] are the Gel'fand-Dikii potentials' 91 . Changing the coupling T k 
induces a flow equation 

du 



dT k 



R k+1 [u], (1.2) 



which is the generalized KdV equation. This equation suggests a possible existance 
of solitons in the flows of these theories. 

Recall that there exist soliton solutions to the KdV equation' 101 

dtu(x, t) + u"'(x, i) — 6u(x, t)u'(x, i) = 0. (1.3) 

Moreover, this equation also has multiple soliton solutions. For the multiple so- 
lutions, the solitons are able to pass through one another and are not necessarily 
of the same size nor moving at the same velocity. Likewise, the modified KdV 
equation 

d t u(x, t) + u'"(x, t) - 6u 2 {x, t)u'(x, t) = (1.4) 

has solutions with these same properties. 

But these solitons are disturbances that propagate over the background, u(x) = 
0. The flow equations that arise in matrix models certainly do not have a limit 
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u(x) = as Tfc — > ±00. For instance, the authors of ref. [11] investigated the 
flow away from the m = 3 multi-critical point by turning on the potential T2. In 
the limit T2 — > — 00 they found that the theory flows to the m = 1 critical point, 
while for the other sign, the theory becomes highly oscillatory. Neither limit has 
u(x) = 0. It might be possible to see a solitary wave propagating over a nontrivial 
background, but this does not occur for the case described in [11]. Instead, as 
Tfc — > 00, infinitely many oscillations develop in a dispersive wave. 

In this paper we will describe a matrix-model which has flows with finitely 
many solitary waves. The model that we will study is a unitary matrix model with 
a potential tuned to its lowest multicritical point, which is essentially QCD on a 
single plaquette [13 ' 141 . To this model we will couple M flavors of nonpropagating 
bosonic quark terms. Integrating out the quarks will lead to an effective action 
with logarithmic terms. The effective potential will have some dependence on the 
rescaled quark mass term, /i. We will show that the square root of the specific 
heat, f(x), will satisfy a fourth order differential equation, where x is the rescaled 
coupling. This equation is related to the mKdV analog of the Gel'fand-Dikii equa- 
tion, which generates the KdV hierarchy. Moreover, we will also show that f(x) 
satisfies a flow equation when \x is varied. We will find that the flow equation is 
not quite the mKdV equation, but the mKdV equation with a dissipative term [lj] 
that is not galilean invariant. 

By studying these equations numerically, we will see that the flow equation 
does indeed describe the propagation of soliton like objects. We have found explicit 
numerical solutions for one and two flavors of quarks, and the corresponding flow 
equations describe propagation of one and two solitary waves respectively. This 
will lead us to conjecture that for M flavors of quarks ther eill be M solitary waves 
in the flow. Unlike the typical soliton solution that propagates over a background 
of u(x,t) = 0, the waves we find propagate over a background that is a solution 
to the Painleve II equation. The solitary waves die out as \i = 00 and move off to 
x = 00 when // = 0; they leave a nontrivial solution to the Painleve II equation 
in their wakes. This will lead us to conclude that the massless quark case has 
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the same specific heat as the infinitely massive case. This remarkable "duality" is 
nonperturbative in nature, that is to say, it is not evident from the equations of 
motion. Note that from hereon, we will refer to these solitary waves as solitons, 
even though they don't quite meet the Zabusky-Kruskal criteria for such objects. 

We have also studied the limit M — > oo to see whether or not it can serve 
as a reasonable regularized version of pure 2d gravity' 16 ' 181 . We will show how to 
calculate the nonperturbative contributions to the specific heat and will show that 
they swamp the perturbative behavior found in the relevant scaling regime. 

In section 2 we derive the differential equation for f(x), as well as its flow 
equation. In section 3 we present our numerical results. In section 4 we study the 
nonperturbative contributions to the specific heat. In this section we also show 
that 2d gravity cannot be reached from a large M theory. In section 5 we present 
our conclusions. 



2. Derivation of the Equations 

In this section we derive a differential equation for f(x), the square root of the 
specific heat, where x is the rescaled coupling. The derivation uses the method of 
orthogonal polynomials, using techniques described by Periwal and Shevitz' 61 and 

[12] 

by Neuberger 

Let us assume that we have a unitary matrix model with M flavors of bosonic 
"quarks" coupled to the N x N matrices. The quark terms are N dimensional 
vectors, ipi and Xii that transform canonically under unitary transformations. Let 
us choose the potential term to be of the form 

where U is a unitary matrix. Since the potential term is quadratic in ip and x, 
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these fields can be integrated out to give the effective potential 

Mtrlog( 7 + I/7 -[/-[/+), 



where we have dropped unimportant constant terms. Notice that potential be- 
comes unbounded below as 7 — > 1 1 1 , hence we expect critical behavior to occur for 
7 near these values. The complete potential for the lowest multicritical point is 
now given as 

V = Ngtr(U + C/t) + Mtr log( 7 + I/7 + U + U ] ). (2.1) 
The partition function is 

Z = J DUe- y V+ u ^ (2.2) 

and is reducible to the diagonal form 

Z= I ^\\\zi- z 3 \\~^ v ^l^ . (2.3) 
J Zmzi . . 

Following [6] and [12], we define the polynomials 

n-l 

P n (z)=z n + Y,a k ,nz k (2.4) 

where all aj-^n are rea l an d the P n 's satisfy an orthogonality relation 

j> dfiP n (z)P m (l/z) = h n 5 n , m , (2.5) 

dz 

d ^ = 2m~z e *v(~ N 9 V ( z + l l z ))- ( 2 - 6 ) 
P n satisfies the recursion relations 

P n+1 (z) = zP n (z) + R n z n P n (l/z) (2.7) 

1 „ , , Z „ , , i?n-2 n+1 



, P n+ l(z) = —P n (z) + ^ Pn+l(l/z) (2.8) 

where ao, n = Rn-i and h n+ i/h n = 1 — i?^. Let us define A n = a n -i jn and 
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B n = a n ,n+2- The recursion relations (2.7) and (2.8) lead to the relations 

A n = A n ^i + R n R n -i, (2.9) 

B n = 5 n -l + i^+l^-l + RnAn-i), (2.10) 

which will result in simple differential equations when we take the limit N — > oo. 
The next step is to consider the integral 



j d^e- v ^ +l '^d{z a P n+l {z)P n {l/z)) 
for the cases a = 0, 1, 2. Integration by parts leads to the equations 

(n + l)(h n+1 -hn) = -j dptV'(z + l/z)(l - l/z 2 )P n+1 (z)P n (l/z), (2.11) 



A n h n = - jdfiV'{z + l/z){l - l/z 2 )zP n+1 (z)P n (l/z), (2.12) 



(2.13) 



((n + 2)(R n+1 R n _ 1 (l - R 2 n ) -A 2 n + 2B n _^j h n = 

- jdnV'{z + l/z){\ - l/z 2 )z 2 P n+1 (z)P n (l/z). 

Now in our case, we find that 

M 

V'(z+l/z) = Ng + — . (2.14) 

v 11 7 + I/7 + 2 + I/2; 

It is clear that in order to derive a differential equation from (2.11)-(2.13) we will 
need to arrange the terms so that we cancel out the z + 1/z term that appears in 
the denominator. This is done by multiplying (2.12) by 7 + I/7 and adding it to 
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the sum of (2.11) and (2.13). After some straightforward but tedious algebra we 
are left with the equation 



2£ n _i - A 2 n + (7 + 1/7) A, + in + 2)(R n+1 R n - 1 (l - R 2 n ) - R 2 n ) + R 



Ng 



—R n (R n+ l + Rn-i) + R n+ 2R n -l + -Rn-2-Rn+l 



— i? n _i_R n+ i (i? n _2-Rr I -l + Rn-lRn + R-nRn+1 + Rn+lRn+2) (2-15) 



+ (7 + 1/7X1+^+1^-1) 

M(l + J R n _ 1 i? n+1 )(l- J R2). 



We now take the limit N — > oo and set g = 1/2. Consider the rescaled variables 



x = (N — n — M)iV~ 1/3 
R n = f{x)N- 1 '* 



7 = -1 - fiN~ 1/3 
A n = N/2 + g(x)N-^ 3 



2B n -! = A 2 n + 2A n - R 2 n - N - M + 2h(x)N-\ 
After substituting these variables into (2.15), we are led to the equation 
2h - fg = x(ff" - {f) 2 - f) + \ff"" - ff" + \{f"? 

-5/ 3 r+2/ 6 -^(/r-(/ / ) 2 -/ 4 ). 

The rescaled versions of (2.9) and (2.10) are respectively 



g' = -f 2 , 



and 



(2.16) 



(2.17) 



h' = lf-ff" + \U'?- 

The last equation is found by observing that 

~ ^n— 1 = iA-n + Ai— l)RnRn— 1 



(2.18) 



— 2A n -2RnRn-l + {RnRn-l) 2 + 2R n R n -iR n ^\R n ^2- 



Taking derivatives on both sides of (2.16) and making the substitutions given 



by (2.17) and (2.18), one finds the equation 

p^^)' = 4//7<, (2.19) 

where H is defined by 

H — f" — 2/ 3 + 2xf. (2.20) 

H is the Painlevell operator, and it is relatively straightforward to show that for 
potentials leading to higher multicritical behavior, the generic form of H is 

i 

where the Tj are essentially couplings and the IZi are the unitary analogs of the 
Gel'fand-Dikii potentials. Moreover, (2.19) is the generator of these potentials. 

M does not explicitly appear in (2.19), but rather arises as an integration 
constant. Multiplying both sides of (2.19) by (H" — fj?H)/ f and then integrating 
gives the equation 

(H" - /j 2 H) 2 = Af 2 ((H'f - /j 2 H 2 + C). (2.21) 

Notice that when fi is zero, (2.21) has the solution H = C, the generic Painleve 
II equation. The integration constant C can be determined by looking at the 
perturbative solutions to (2.21). In particular, the asymptotic expansion derived 
from (2.21) for large x is 



r i/2 r 1 / 2 ,, 2 



but it is known that for /i = 0, the leading order expansion is 116 181 

f( x ) = x 1 ' 2 - — + 0{x^l 2 ). (2.22) 
2x 

Hence we see that C = Afj, 2 M 2 . 
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We conclude this section by deriving a flow equation for f(x) as a function 
of ii. The flows are basically determined by ward identities of the matrix model 
partition function. We can take a derivative with respect to 7 on both sides of 
(2.5), giving the equation 

-— / dp *\ + Z ^ l [\ . P n (z)P n (l/z) = d^K. (2.23) 
7 J 7 + I/7 + z + 1/z 

We can also do the same thing after inserting z and 1/z into (2.5), giving the 
equation 

Mr 2-f + z + l/z n . . „ . . . 

<P da '—— zPJz)PJl z) 

7/^7 + 1/7 + 2 + 1/2 nv ' ny 1 ' (2.24) 

+ h n d^A n h n R n d^R n ~\ = d^(h n R n R n — 1), 

and an analogous equation for the 1/z insertion. Combinining all three equations, 
we can deduce the expression 



• — (7 - RnRn-i) = (7 + 1/7 - 2R n R n ^)^ - 2R n ^ 1 d 1 R n - 2d T 

7 



Since h n is given by 

n-1 

/i n =^n( 1 -^ 2 )' 

i=0 

we can derive the recursion relation for its derivative 

9 7 /l n+ i - — d^hn = — ^— ( 2 - 26 ) 

Taking a difference of (2.25) and invoking (2.9) and (2.26) leads to the equation 
2M (R n+1 - i? n _i) =(7 + 1/7 - 2R n+1 R n ) ~ 2R ^ R - 



— 2Rn(Rn+i — Rn-i)—r 2i? n <9 7 (_R n +i + R n -\). 

(2.27) 

We can get rid of this last d^h n /h n term by dividing (2.27) by R n (R n+ i — R n -i), 
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taking another difference and then multiplying by R n +i — R n -i- Hence we find 
2R n (R n+ i — R n ^i)d 1 R n R n +i — Rn-i 







1 — Rn Rn+2 — Rn 

x ( (7 + 1/7 - 2R n+2 R n +i)p^ ± - + d^Rn+2 + Rn)) (2.28) 

V 1 - R n +i / 



L n+1 

+ (7 + I/7 - 2R n+1 R n )^^ + ^(^n+l + Rn-l) 



'"n 



Taking the double scaling limit then leads to the equation 

= Aff'drf + fi%f + %ff'd u f+Af 2 d a f - d u f" 



which, after dividing through by /' and changing the order of derivatives, gives 



Hence the flow equation is defined by a linear operator and the flows are determined 
by finding the eigenfunctions with zero eigenvalue. 

Now let d u f = ip. If we substitute J ip for H in (2.19), then it is straightforward 
to show that the equation reduces to the flow equation, (2.30). Therefore, we find 
that 

duf = CH', (2.31) 

where H is given by (2.20). The constant C can be determined by looking at the 
asymptotic expansion for /. The leading behavior is easily determined from (2.21), 
giving 

f _ q M | M/i 2 | 
2x 16x 2 

and hence the leading term for d n f is 2M/x/16x 2 . The leading term for H' is 
Mfj, 2 /Ax 2 , thus C = l/2/i. Letting t = log//, we can then write the flow equation 
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as 



(2.32) 



Equation (2.32) is almost the modified KdV equation, the difference being 
that there is extra term, xf + /. This term does two things — it breaks galilean 
invariance and it dissipates wave motion. Nevertheless, this equation has soliton 
like solutions, albeit solitons that die out. Our numerical results will exhibit this 
behavior. 

3. Numerical Results 

A relaxation method was used to find solutions to (2.21) with the asymptotic 
behavior given by (2.22). Following the work of [19,11] we consider the diffusion 
equation 

d T f = D(f), (3.1) 

where D(f) is 

£>(/) = {H" - fi 2 H) 2 - 4/ 2 ((#') 2 " V 2 (H 2 - AM 2 )). 

If / relaxes to a r independent solution, then we have found a solution to (2.21). 
One starts with a test function with the desired asymptotic behavior and lets it 
evolve. This is done by using the recursion 

/new = /old + hD(f i d ). 

However, h cannot be chosen too large, otherwise the system destabilizes. 

The problem with this approach is that a solution is not always reached, some- 
times / blows up. In fact for /i 2 less than a certain value it turned out to be 
nearly impossible to find stable solutions. However if we can find a solution for a 
given value of fj 2 , we can then employ the flow equation (2.32) to find approximate 
solutions for smaller ji 2 . 
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A starting function was generated by considering solutions to the Painleve II 
equation 

g " = 2g 3 - 2xg + 2M. (3.2) 

It turns out that for M < 1/2 there is a solution to (3.2) with asymptotic behavior 
g ps H-y^ as x — > oo, which is free of poles. The asymptotic behavior of this solution 
as x — > — oo is 5 « M/x. M was set to M = 0.49 and a numerical solution was 
found by using the relaxation method, where the starting function was chosen to 
be a solution to the algebraic limit of (3.2). The algebraic solution with the desired 
asymptotic behavior has a jump at the critical point x = 3(M/2) 2 / 3 , but this did 
not prove to be a hindrance in finding a solution to (3.2). 

The resulting solution of (3.2) was then used as a starting function for (3.1). 
Since the asymptotic behavior of a solution to (3.2) is slightly different than that for 
(2.21), the endpoints of f(x) were allowed to move for a large number of iterations, 
but were eventually fixed at x — ±20.0. The difference in asymptotic behavior is 
most apparent in the strong coupling regime (x — > — oo) where the behavior for 
nonzero /i is / ~ e^ x . With M still set to 0.49, we were able to reach a solution 
for (j, — 1. This solution was then used as the starting function for different values 
of M and /i. 

Figure 1 shows the results for M = 1. The solution for /i — 2.1 was found using 
relaxation and the lower values of /i were reached using the flow equation. The 
flow was calculated using a fourth order Runge-Kutta routine. Moreover, after 
each step in the flow, a conjugate gradient routine 1201 was used to minimize the 
value of Y^i{D(f(xi))) 2 . We used two different lattice spacings (Ax = .3 and .1) 
and two different flow steps (At = —.01 and —.001) with no qualitative difference 
in results. The flow is glacial because of the log dependence of the values of 
li plotted here start at fj, — 2.1 and end at /i = 1.1 x 10~ 3 Figure 1 also has a 
plot of a solution to (3.2) with M = and weak coupling behavior —yfx. Since 
the specific heat of the theory is given by — / 2 (x), it is clear that as fi approaches 
zero, the behavior of f(x) in the strong coupling regime approaches that of pure 
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QCD. The only region where f(x) differs markedly from this behavior is in the 
region where f(x) increases sharply. Let us contrast this with the behavior of f(x) 
as /x — > oo. In this case the quark terms should decouple and the theory should 
also behave like pure QCD. In this case f(x) approaches the solution of (3.2) with 
M = but with asymptotic behavior y/x. Hence we see that specific heat in the 
scaling regime is identical for the zero-mass and the infinite mass case. 

Let us next compare the specific heat for the case M = 1 and assorted values 
of /i with the specific heat of pure large iV QCD. Figure 2 shows the difference in 
specific heats, /qqd (x) — f 2 (x) , for values of /x ranging from 2.1 down to 1.1 x 10~ 3 . 
The figure suggests a soliton propagating in the negative x direction as the value 
of /x increases, but with dissipating size. When /x — 0, the soliton is off at positive 
infinity with no wake in the finite x regime. As the wave travels in the — x direction 
it eventually dies out, leaving pure QCD behavior for the specific heat. 

This dual structure between the two mass extremes is quite surprising, at least 
from the standpoint of the original matrix model partition function. It is a non- 
perturbative effect and is not seen in a saddle point calculation of the eigenvalues. 
The peak that appears in the specific heat is similar to the Schottky anomaly, 
which is the result of tunneling in a two well system. If the first N — 1 eigenvalues 
relax to the ground state, then the N th eigenvalue will have an effective potential 
with two wells, provided that /x is small enough. The presence of this second well 
will result in a relatively large density of states around some corresponding energy, 
leading to a peak in the specific heat. The effects of these wells will be discussed 
further in the next section. 

We have also succeeded in finding a solution for M = 2 and /x = 3.5, although 
only for lattice sizes larger than or equal to Ax = 0.3. With this solution we then 
flowed to lower values of /x, giving the results shown in figure 3. Here we see that 
an extra bump develops as /x gets small. Unfortunately, we were unable to reach 
values of /x below .01 without developing large errors in the value of J2i(D(fi)) 2 . 
However, this seemed far enough to clearly see the second peak. The peaks are 
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quite clear in figure 4, where we have graphed /qcd — for values of \x between 
3.5 and 0.1. In this case there appear to be two soliton peaks where one clearly 
propagates in the — x direction and eventually dies out. The second peak definitely 
moves slower than the first, although it is not clear from our data if it even moves 
at all. However, we deem it likely that this peak eventually moves off to infinity 
as /i approaches zero. 

It is natural to postulate that a system with M bosonic flavors will have M of 
these peaks in its specific heat. It is reasonable to expect that the peaks are due to 
classical configurations corresponding to some sort of instanton processes, which 
are clearly a result of the quark terms in the effective action. A question worth 
asking is whether these configurations are also present in two dimensional lattice 
QCD with propagating quarks. Recall that the theory we have been considering 
is QCD restricted to one plaquette. If there are no quarks in the theory, then 
this is equivalent to QCD on a lattice, since in the full lattice theory the plaqettes 
decouple. But once propagating quarks are introduced the plaquettes become 
coupled. However it is possible that the peak structure found here will have some 
remnant in the full two dimensional theory, if, for example, the two dimensional 
theory has soliton solutions that are localized in space. This problem is left for 
future work. 

4. Gravity is Unreachable 

In [16] it was argued that as the number of flavors M becomes large, then 
under a rescaling, the Painleve II equation reduces to the Painleve I equation. This 
suggests that pure two dimensional gravity can be reached from two dimensional 
QCD. However, in [18] it was shown that the theory with the correct asymptotic 
behavior will have poles in f(x) if we insist that f(x) is real. This was shown 
to contradict the Schwinger- Dyson equations. It was then suggested that under 
a suitable regularization we could remove the poles and reach pure gravity. The 
poles were arising because the effective matrix potential was unbounded below. 



14 



This is a consequence of setting p to zero. If we were to bound this by introducing 
a nonzero p, then we can consider the simultaneous limit p — > and M — > oo and 
check whether gravity is reached. We will now show that this won't occur. 

When M is large we can set e = 2/M and define x and f(x) such that 

x = 3e -2/3 + ( 3 /4) l/5 e 2/15 5 

(4.1) 

/ = e -V3 + (1/18) l/5 e l/15 / ~ 

Then to leading order in e the Painleve II equation reduces to Painleve I. Therefore, 
in order that the regularized theory reach pure gravity, the solutions that satisfy 
the boundary conditions must pass sufficiently close to the critical point x = 3e~ 2//3 
and / = e -1 / 3 . 

Let us next consider the nonperturbative behavior that can be derived from 
(2.21). Define f(x) as fasy(x) + g{x), where / a sy is the asymptotic expansion in 
1/x and g(x) is a subdominant term. To leading order in its asymptotic expansion, 
H(x) is given by 

H(x) = 2M- Mp 2 /Ax + ... + g"(x) - Axg(x) + ... 

(4.2) 

= 2M - M/i 2 /4x + ... + h{x) 
Plugging this into (2.21) leads to a linear equation, which to first approximation is 

-4/j?Mh" = 4(x - M/x 1 / 2 )((M/2 2 /2x)h' - 4M/i 2 /i). (4.3) 

We have assumed that \i is small and thus higher order terms in \i can be dropped. 
Letting h = p{x)e~^/^ x 1 , we can use a WKB approximation to derive the equa- 
tion 

4p' + p = -p + 4Mp, (4.4) 
thus giving p = Cx^ 4 " 1 / 2 , where C is a constant in x that depends on M and p. 
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Now g(x) satisfies the approximate equation 

h(x) = g"{x) - (Ax - (6M/x^ 2 ))g(x). (4.5) 
This leads to g(x) = r(x)e~^ 3 ^ x ' where r(x) approximately satisfies the equation 
p(x) = -Ax l / 2 r'{x) - r(x)/x 1/2 + QMr(x)/x 1 / 2 . (4.6) 

Hence r(x) can be approximated as r(x) « Cx M /(2M — 1) 

We can derive the /i dependence of C using the flow equation. From (2.32) we 
find that 

dtg = \ti w-(2M-l)y, 
from which we derive that C ~ e -( 2M -^ = l/pi 2M -\ and 

^ ^(^ 2 ) M e -(4/3),3/ l 

yv ; 2M-1 v y 

k depends only on M. From this form of g(x) we see that when /j, — > 0, the 
perturbative expansion will be wiped out by nonperturbative effects, if M > 1/2. 
However, for finite M, if \x is small then there is a large range of values for x where 
the perturbative expansion of the Painleve II equation is accurate to large order 
in 1/x. This is because the nonperturbative piece only depends on a power of //. 
Hence, while nonpertubative effects ruin the limit to Painleve II, the effect is quite 
gentle. 

If we now consider the large M limit, it is possible that these nonperturbative 
effects will be gentle enough that they won't destroy the Painleve I behavior, 
perhaps because the overall coefficient in front of the perturbative piece is small. 
That is to say, we can choose \i to be small and yet have f(x) pass sufficiently close 
to the critical point. We now demonstrate that nonperturbative effects do indeed 
destroy the limit. 
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In order to understand this, let us consider the origin of the nonperturbative 
contribution for the lowest multicritical point of the unitary matrix models. The 
square root of the specific heat is a solution to the Painleve II equation with M — 0. 
The perturbative behavior can be determined from the distribution of eigenvalues 
on a circle. The potential is given by 

V{9) = 2Ng cos0i - ]T log(sin 2 \{9, - 9j)), (4.8) 

i i<j 

and the equations of motion are 

= -2gNsm6i -iV^cot±(^ -0j). (4.9) 

Defining x = i/N,dx = l/N gives the equation 

l 

= -2gsin9(x) - P J dxcot \{9{x) - 9(x')), (4.10) 
o 

where the P in front of the integral indicates the principle value. Introducing the 
density of eigenvalues u(X) = dx/dX, leads to the equation 

a 

= -2gsmd - J d\cot\(d - A), (4.11) 

—a 

where ±a are the end points of the eigenvalue distribution. 

The density of states can be determined in a way analogous to the method 
used in [21]. We define the function F{9), 

a 

F(9) = J dX cot \(9 - A), (4.12) 

—a 

which is analytic everywhere on the strip — n < Re9 < n, except for a cut on 
the real line between —a and a. Using the equations of motion and the fact that 
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F(6) — > =i=i as Im# — > ±00, results in 

U (A) = -a/2(cosA - cos a) cos(A/2). (4.13) 

71 

a is related to g by 

#(l-cosa) = l. (4.14) 

This density of eigenvalues then gives the perturbative expansion to the theory. 

However, it is clear that the potential has a saddle point if one of the eigenvalues 
is sitting at 6 = ir. Hence we look for a solution where an eigenvalue is fixed at 
this position and the others are allowed to seek a local minimum. The potential 
for this configuration is given by 

V(6) =2gN cos 6 -^log(sin 2 \{6 - 6$) 

^ (4.15) 
+ 2NgJ2^s6i- J2 log(sin 2 i(^-^)), 

i^l Ki<j 

where 6 = ir. But this is the potential for a U(N— 1) system with one fermion flavor 
and we know what its perturbative expansion looks like. Hence we can determine 
the nonperturbative contribution by comparing the perturbative behavior of these 
two theories. 

The partition function for a unitary matrix model is given by 

N 

Z N = e- FN = N\ Y[{1 - Rf) N+1 -\ (4.16) 
i=i 

If we assume that R? = i/N, which is the tree level contribution, then F/v is given 
by 

F N = - \og(M) - - i) log(l - i/N) 

1 o (4.17) 

« - log AH - 1 ; (log - 3/4) + C. 
SV ; 2gN K &V 2gN ' 1 > 

Comparing this with F/v-i, we find that these two free energies differ from each 

other by log(iV) — 1 if f? is in the double scaling region. But we also need to take into 
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account the difference in the free energy from the scaling region. This contribution 
satisfies F" = f 2 (x). For the zero quark case, the perturbative expansion of f(x) 
is 

f(x)=x 1 / 2 -^ + ..., (4.18) 
while for the one fermion case, f(x) is 



m = ^ + Mil - *"V« : + 1/16 + 



with M — 1. These expansions are each easily found from the Painleve II equation. 
The difference in f 2 (x) for the two situations is then 



6f\z) = JL-W + ... (4.20) 



and therefore the difference in the free energy is 

AF S = (4/3)x 3 / 2 + \ Xo^xN- 2 ^ + ... (4.21) 

The appearance of N^ 2 ^ is the usual situation of unsealed variables appearing in 
the logarithm terms. 

Given the above, we see that the nonperturbative contribution to the partition 
function is the perturbative contribution multiplied by the factor 

c ^7^^ (4/3,i3/ 7*" H '" w/l (4 - 22) 

where C is a constant of order unity. The log iV term that appears in the difference 
of the two free energies disappears because any one of the iV eigenvalues can sit 
at the top of the potential. Strictly speaking, the density of eigenvalues used in 



19 



evaluating V" should be that for the M — — 1 system. For a nonzero value of M, 
u(X) is modified to 

If (M/7V)( 7 i/2 + i/ 7 V 



*A (7 + 1/7 + 2 cos A) ^7 + I/0 + 2 cos a/ (4.23) 
x a/2(cosA — cosa) cos(A/2). 

However, the contribution coming from the term multiplying M will be suppressed 
by a factor of 1/x 3 / 2 , so we can ignore it. Taking two 6 derivatives of (4.11) gives 



V"(n) = -2 9 N + ^ / ^ ^(cosA -cosa) ^ 

K ' y 2tt J cos(A/2) v ; 



The integral that appears in (4.24) can be evaluated, satisfying 



a 

I 



„ \/2(cos A — cosa) . 

rfA cos(A/2) U 'ft -«"(»/»)). <«5) 



and hence leading to V"(ir) = 2gNcos(a/2). Letting g = (1 + a;iV _2 / 3 )/2 and 
using (4.14), we find that V"(7r) = -x l / 2 N 2 / z . Since V" is negative, the integral 
in (4.22) is not bounded below. But in any case, the limits of integration should be 
bounded such that the argument of the exponent is not much greater than unity, 
since the eigenvalue will eventually meet a repulsion from the other eigenvalues. 
Hence the nonperturbative contribution to the free energy is 



Fnp = Cx-^e-W^ , (4.26) 
and thus its contribution to f(x) to leading order is 

f np = Cx-^e~W^ 3, \ (4.27) 
This is precisely the behavior that one would find for a WKB solution of (3.2). 
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As an aside, we note that for the hermitian case, if there is a saddle point 
such that one eigenvalue is sitting at the critical point, then the same analysis 
could be applied to determine the nonperturbative behavior. In particular, this 
would imply that the free energy of the disk plays some role in this behavior. For 
the unitary case, the dominant contribution to the nonperturbative piece comes 
from the free energy of the disk, which is proportional to 1/g, where g is the 
renormalized coupling. But actually, the hermitian matrix models either have no 
saddle points, which is the case for the odd potentials, or they are unbounded below. 
However, it is true that for the hermitian case, the strength of the nonperturbative 

[22 23] j [24] 

piece ' is proportional to e~ > g . Moreover, it was shown by David that the 
nonperturbative piece for the pure gravity case is a consequence of an eigenvalue 
sitting at the top of a potential. Hence it is possible that the free energy of the 
disk will play some role in determining the nonperturbative behavior of a generic 
string theory. We won't have anything more to say about this subject. 

Now let us turn to the situation with M flavors of bosons and let us assume 
the \i is small compared to 1/x. If /i is small enough then the saddle point will 
actually be a local minimum. The effective potential for one eigenvalue fixed at 
9 = 7r is 

V(9) =2#iVcos# + log(7+ I/7 + 2 cos + (M — 1) ^log(sin 2 \ {9 — Oj)) 

+ 2NgJ2^s6 l - log(sin 2 i(^-^)), 

i^l Ki<j 

(4.28) 

which is the potential for a theory with M — 1 flavors, if 9 is fixed at n. The differ- 
ence in free energies between these two theories is again found from the differences 
in f 2 (x), which is 

Af 2 (x) = (M - (M - l))/x 1/2 + (M 2 — (M — lf)/2x 2 + ... 
= l/a: 1/2 + (M-l/2)/a; 2 + .... 
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From this we find the difference in free energies to be 



AF S = (4/3)x 3 / 2 — (M — 1/2) log(:riV- 2 / 3 ) + .... (4.30) 

We should also add the term M log(/i 2 iV~ 2//3 ) to AF S , which is the contribution to 
the potential from the eigenvalue sitting at 6 = ir. The main contribution to V"(n) 
comes from the second term on the r.h.s. of (4.28), thus V"(ir) w 2/(/z 2 iV~ 2 / 3 ). 
Unlike the previous example, V" is positive and thus this eigenvalue sits in a well. 
Hence we should expect peaks to appear in the specific heat. Putting this all 
together, we find that the nonperturbative contribution to the free energy is 

~M-l/2 

F - C— e" (4/3)a; (4 31) 

p ~ mvv m+1 

and therefore the contribution to f(x) is 

f - c - e" (4/3)a; (A 32) 

Inp ~ M 1 / 2 /2 2M + 1 ' 

Hence the behavior matches the subdominant contribution to the asymptotic ex- 
pansion of (2.21). 

It is clear from these examples that in order for the nonperturbative behavior 
to avoid swamping the perturbative expansion, it is necessary that AF S diverges 
when the scaled variable, in this case x, becomes large. At the same time we still 
want to remain in the scaling regime. To this end, consider what happens in the 
large M limit, with the rescaling of variables given by (4.1). Given an x there is a 
corresponding x, depending on the value of M. If we fix x, but decrease M by 1, 
then x shifts by x — > 5 + 2(4e/3) 1 / 5 . i.e. as M — > oo, the shift in cc approaches zero. 
Now in order for the limit M — > oo to be able to approach 2d gravity, it needs to 
have a sensible perturbative expansion in x. However, if we assume that shifting 
M by 1 cannot affect this perturbative expansion, and since the shift only changes 
x by an infinitesimal amount, then it must be the case that a nonperturbative 
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contribution coming from an eigenvalue sitting at 9 = tt must contribute a piece of 
order unity if fi ~ 1. In other words, the nonperturbative piece is not suppressed 
by a factor of e 1//ff , where in this case, g = x~ 5 / 4 . Hence it must be true that as 
H — > 0, the curve f(x) must be far from the critical point and therefore 2d gravity 
won't be reached using this regularization. 

5. Discussion 

As we have already mentioned at the end of section 3, an important question 
is to what extent the results presented here extend to the full two dimensional 
lattice QCD. For instance, we might wonder whether these solitons remain in this 
extended theory. If they correspond to localized objects, then it would seem likely 
that there is a remnant in the full 2d theory, because their structure appears at the 
single plaquette level. If there is a remnant in the two dimensional theory, then it 
might turn out that there is an interesting relation between the solitons one finds 
in the flows and the solitons that arise as nontrivial solutions of the classical field 
equations. 

There is also a question of integrability. We have been rather cavalier in re- 
ferring to these solitary waves as solitons. In its strict sense, the name soliton, 
coined by Zabusky and Kruskal, implies that two such objects will pass through 
one another, such that out at infinity, the only signature that they leave on one an- 
other is a phase shift. This is a direct consequence of the integrability of the KdV 
equation. For the waves described here they die out while they are passing through 
each other, so it is hard to see if there are infinitely many conserved charges in 
the theory. However, it is possible that there is a notion of integrability in this 
system and an inverse scattering transform can be found to solve the problem of 
wave interaction. Perhaps there is a change of variables where such a structure 
would be made more transparent. In any event, the fact that this flow equation 
shows up in a solvable model lends credence to this supposition. 
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FIGURE CAPTIONS 



1) Flows of f(x) for M — 1 and different values of fx ranges from 2.1 down 
to 1.1 x 1CT 3 . Also shown is a plot of f(x) for M — 0. 

2) Plots of /q CD (x) - with M = 1 and ranging from 2.1 (bottom) to 
1.1 x l(T 3 (top). 

3) Flows of /(#) for M = 2 and different values of /i. /x ranges from 3.5 down 
to 1(T 2 . 

4) Plots of fq CB (x) - with M = 2 and fi ranging from 3.5 (bottom) to 
1.1 x 10 _1 (top). 
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